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Using two difTerent approaches known in the hterature, both based on self-adjoint extension 
method, we propose a general approach, which is based on the conditions imposed by the physical 
problem. One of the advantages is that this approach yields the self-adjoint extension parameter in 
terms of physics of the problem. We apply it for the spin-1/2 Aharonov-Bohm problem in conical 
space in the nonrelativistic limit. The self-adjoint extension parameter to the bound state and 
scattering scenarios are determined. Our proposal solves the problem of the arbitrariness of the 
self-adjoint extension parameter, proposed in Ref. [Phys. Rev. D 40, 1346 (1989)]. The present 
method is general and suitable for addressing any quantum system with a singular Hamiltonian that 
possesses bound and scattering states. As an application, we apply it for the spin-1/2 Aharonov- 
Bohm problem plus a two-dimensional isotropic harmonic oscillator. 
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I. INTRODUCTION 

The Aharonov-Bohm (AB) effect [1] is one of most 
weird results of quantum phenomena: a charged parti- 
cle feels effects due non accessible magnetic field only 
travelling at a magnetic force-free region. The charged 
particle wave function then acquires a phase shift which 
would yield produce measurable interference effects. This 
occurs because the magnetic potential is not zero in the 
force-free region. Recently it was reported that the AB 
effect is also realized for neutral particles like photons 
[2]. In the AB effect of spin-1/2 particles [3, 4], besides 
the interaction with the magnetic potential, an additional 
two dimensional 6 function appears as the mathematical 
description of the Zeeman interaction between the spin 
and the magnetic flux tube [5, 6]. This interaction is the 
basis of the spin-orbit coupling, which causes a splitting 
on the energy spectrum of atoms depending on the spin 
state. Hagen [7] argued that this 6 function contribution 
to the potential can not be neglected when the system 
has spin, having shown that changes in the amplitude 
and scattering cross section are implied in this case. 

Point interactions usually appear in quantum systems 
in the presence of topological defects and have been of 
great interest in various branches of the physics for its 
relevance as solvable models [8] . A simple but non trivial 
example is the case of a cone arising from an effective 
geometry immersed in several physical systems, such as 
cosmic strings [9-14], defects in elastic media [15], de- 
fects in liquid crystals [16], and so on. In such systems, 
although the particle does not have access to the core 
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(defect) region, its wave function and energy spectrum 
are truly influenced by it. 

In quantum mechanics, singularities and pathological 
potentials are in general dealt with some kind of reg- 
ularization procedure. A common approach to ensure 
that the wave function in the presence of a singularity is 
square-integrable (and therefore might be associated to 
a bound state) is to force it to vanish on the singularity. 
More appropriately, an analysis based on the self-adjoint 
extension method [17] encloses the boundary condition 
possibilities that still give bound states. The physics of 
the problem determines which of these possibilities is the 
right one, leaving no ambiguities [3, 18]. This method 
has been applied by several authors, in particular, for 
AB-like systems [3, 6, 19-22]. However, these works pro- 
vide the most important results (e.g., energy spectrum, 
phase shift, S matrix) in terms of an arbitrary real pa- 
rameter, the so called self-adjoint extension parameter. 

The self-adjoint extension of symmetric operators is 
a very powerful mathematical method end it can be ap- 
plied to various systems in relativistic and nonrelativistic 
quantum mechanics, supersymmetric quantum mechan- 
ics and vortex-like models. The AB interaction within 
the framework of k deformed Pauli equation was studied 
in Ref. [23]. The structure of additional electromag- 
netic flelds to the AB fleld, for which the Schrodinger, 
Klein-Gordon and Dirac equations can be solved exactly 
was addressed in [24]. Correa et al. [25] analyzed the 
supersymmetrics of a spin-1/2 particle in the context of 
magnetic vortex and anyons. They found that there are 
just two self-adjoint extensions for the Hamiltonian that 
are compatible with the standard N = 2 supersymme- 
try. Other approaches for the AB problem supersymmet- 
ric extension and for the origin of the hidden supersym- 
metry were addressed in Refs.[26] and [27], respectively. 
Slobodeniuk et al. [28] presented a comparative study 
of the electronic properties of a two-dimensional elec- 
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tron gas in the presence an infinitesimally thin magnetic 
field solenoid with relativistic dispersion in graphene and 
quadratic dispersion as in semiconducting heterostruc- 
tures. The same study was performed for a graphene 
system in the inhomogeneous magnetic field consisting of 
the Aharonov-Bohm flux and a constant background field 
[29]. In Rcf. [30] the Aharonov-Bohm-Cashcr problem 
was examined for a charged particle describing a circular 
path in the presence of a Lorentz- violating background 
nonminimally coupled to the spinor and gauge fields. For 
other systems where it is possible to apply the regular- 
ization method, see [31-33]. 

In this paper, we present the details of a general physi- 
cal regularization method recently proposed in Ref. [34] . 
It is based on the physics of the problem and one of his 
particularities is that it gives us the self-adjoint exten- 
sion parameter for both bound and scattering scenarios. 
The description of the formalism is based on the works 
of Kay-Studcr (KS) [35] and Bulla-Gesztezy (BG) [36], 
both using the self-adjoint extension method. We also 
apply the method for the spin- 1/2 AB problem plus a 
two-dimensional isotropic harmonic oscillator. 

The paper is organized as follows. In Sec. II we write 
the Hamiltonian of the spin- 1/2 AB problem and derive 
the equation of motion that governs the dynamics of the 
particle. In Sec. Ill we present the KS and BG self- 
adjoint extension methods used in the formulation of the 
regularization method proposed here. The KS method 
has the advantage of yielding the self-adjoint extension 
parameter in terms of the physics of the problem, but it 
is not appropriate for dealing with scattering problems; 
on the other hand, the BG method is suitable to address 
both bound and scattering scenarios, with the disadvan- 
tage of allowing arbitrary self-adjoint extension parame- 
ters. Further, we also derive the expressions for the en- 
ergy bound state, phase shift and the scattering matrix 
in terms of the physics of the problem. By combining the 
KS and BG methods, a relation between the self-adjoint 
extension parameter and the physical parameters of the 
problem is found. In Sec. IV we apply the method for 
the spin-1/2 AB problem in a two-dimensional isotropic 
harmonic oscillator. We derive the expression for the par- 
ticle energy spectrum and analyze it in the limit case of 
the vanishing harmonic oscillator potential recasting the 
result of pure spin-1/2 AB problem in conical space. In 
Sec. V we present a brief conclusion. 



II. THE EQUATION OF MOTION 

The idealized situation of a relativistic quantum par- 
ticle in the presence of a cosmic string is an example of 
a gravitational effect of topological origin, where a parti- 
cle is transported along a closed curve around the cosmic 
string [37]. This situation corresponds to the gravita- 
tional analogue of the electromagnetic AB effect with the 
cosmic string replacing the flux tube [38-42] . Such effects 
are of purely topological origin rather than local. The 



bound state for the spinless AB effect around a cosmic 
string was addressed in [43]. The authors observed that 
the self-adjoint extension of the Hamiltonian of a par- 
ticle moving around a shielded cosmic string gives rise 
to a gravitational analogue of the bound state AB ef- 
fect. Here, our initial proposal is to analyze the spin-1/2 
AB problem in the cosmic string spacetime with an in- 
ternal magnetic field. The cosmic string background is 
described by the following metric in cylindrical coordi- 
nates {t, r, ip, z): 



-dt^ + dr^ + a^r'^dip^ + dz^ 



(1) 



with -cxo < (t,z) < oo, r > and < < 27r. The 
parameter a is related to the linear mass density to of 
the string by a = 1 — 4to runs in the interval (0, 1] and 
corresponds to a deficit angle 7 = 27r(l — a). The ex- 
ternal gravitational field due to a cosmic string may be 
approximately described by a commonly called conical 
geometry. Usually, only the case a < 1 is considered in 
cosmology, since a > 1 corresponds to a negative mass 
density cosmic string. For a = 1, the cone turns into 
a plane. The spacetime produced by a static straight 
cosmic string can be obtained in the weak-field and low 
speeds limit (valid for m ^ 1) [37]. The above metric has 
a cone-like singularity at r = and the curvature tensor 
of this metric, considered as a distribution, is given by 



R\l = R\^r1^2t: 



1 



(2) 



where 5^{r) is the two-dimensional delta function in flat 
space [44] . This implies a two-dimensional conical singu- 
larity symmetrical in the 2-axis, which characterizes it as 
a linear defect. 

In order to study the dynamics of the particle in a non- 
flat spacetime, we should include the spin connection in 
the differential operator and define the respective Dirac 
matrices in this manifold. The modified Dirac equation 
in the curved space reads (ft = c = 1) 



(3) 



where q is the charge, M is mass of the particle, '^{x) is 
a four-component spinorial wave function, and is the 
spin connection, which is given by 



r,„ = --i^^^h^'^el 



T 'T '6(^)6(6)^^, 



(4) 



and 7^ — 6(0) (2^)7'"' E^re the 7 matrices in the curved 
spacetime. We take the basis tetrad [45-47], 
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(5) 



(with a = 1 giving the flat space-time) satisfying the 
condition 



(6) 
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with g^'^ = diag(— , +,+,+). For this conical spacetime 
the spin connection can be expressed by 



(7) 



where k = (1 — a)/2ar and 



Y = cosipj^^^ +sin(/57^^' ~ cr'' 'o ) ■ 
Moreover the a matrices are now written as 









( 


1) 




J 









(9) 



where cr' = (ct'', cr"^, ct^) are the Pauli matrices in cyhn- 
drical coordinates obtained from the basis tetrad (5). 

For the specific tetrad basis used here, the spin con- 
nection is 

r^ = (o,o,r^,o), (10) 

with the non- vanishing element given as 



with 
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(11) 



(12) 



We are interested in the non-relativistic limit of the 
Dirac equation, so it is convenient to express it in terms 
of a Hamiltonian formalism: 



(13) 



with 



H = a' (-^9, + qA,) + qAo - ^7"7^^^ + /3M, (14) 
where we have defined the notation 

«' = 7V, 7° = /3, 

and the Dirac matrices obey the anticommutator rela- 
tions 

In the non-relativistic limit we have to extract the tem- 
poral behavior of wave function. Writing 



= e 



-iMt 



(15) 



we obtain the equations 

{a'n' +iKa'') X + qAo'i> = idt<^, (16) 

(crV + ina'') $ + qAoX - 2MX = idtX, (17) 



where tt* = —idi + qAi. In the non-relativistic limit the 
"small" component of the wave function, X, can be re- 
lated with its "large" one, $, by 



^ = ^('^V+Z«:CT'-)$. 



(18) 



Now substituting this result in Eq. (16), we find 

Hnr'^ = tdt^, (19) 

where 



Hnr = — (crV' -t- iKCr'") (cr^TT^ + iKCr'') + qAo- (20) 

The magnetic flux tube in the background space de- 
scribed by the metric above considered is related to the 
magnetic field by 

.B = .(VxA) = f?^z, (21) 
a r 

where (j) = 0/27r is the flux parameter, the vector poten- 
tial in the Coulomb gauge is 



ar 



Ao = 0, 



(22) 



and s = ±1 is twice the spin projection parameter [7]. 
The parameter s implies that the Dirac equation de- 
scribes the planar motion (in the absence of the z co- 
ordinate) of the particle having only one projection of 
the three-dimensional spin vector. The choice (21) also 
gives a flux tube coinciding with the cosmic string and 
the z-axis. 

Using the relations 



_i_7 ii I ■ ijk k 

a a-' = g ■' + le ■' a , 



1 



-5'' 



the Hamiltonian (20) can be written as 
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where 



and 



2ar 



a 



11 5{r) 
2M r 



V 
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" dr"^ r dr a^r^ dip^ ' 



(23) 
(24) 



, (25) 



(26) 



(27) 



is the Laplacian operator in the conical space. For this 
system the total angular momentum operator. 



J = -id^ + -cr^, 



(28) 
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Scattering 
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-1 
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Scattering 
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< 


-1 
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Bound and Scattering 
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> 


+1 


< 


< 


Bound and Scattering 



TABLE I. Possibilities to occur Uai 



> and Us 



< 0. 



commutes with the effective Hamiltonian. So it is possi- 
ble to express the eigenfunctions of the two dimensional 
Hamiltonian in terms of the eigenfunctions of J. These 
functions have the form 

-^£t ( /£(r)e'("-^/2)'^ " 



(29) 



with m = n + 1/2, n € Z. At the same time, the radial 
equation for /g (r) becomes 



where 



H — Ho + ^sho 



^0 = -77V-r 



1 

21/ 



short 



1 d 

r dr 

T] 5{r) 
2M r ' 



and 



J 



-(m - 



^ 1 



-)• 



(30) 

(31) 

(32) 

(33) 
(34) 



The Hamiltonian in Eq. (31) governs the quantum dy- 
namics of a spin- 1/2 charged particle in the conical space- 
time, with a magnetic field B along the z-axis, i.e., a 
spin-1/2 AB problem in the conical space. Let us con- 
sider a conical defect with a radius ro nucleus, so it is 
suitable to write Z^short('') as [5, 7] 



Z^short(r) = 77 



S{r - ro) 



(35) 



and, at the end, the limit ro — >■ is taken. Although 
the functional structure of Z4hort and Z^short are quite 
different, as discussed in [7], we are free to use any form 
of potential provided that only the contribution of the 
form (33) is excluded. 

Let us now study the Z^short signal, which provides 
us the possibilities to obtain both scattering or bound 
states. Therefore, it is important to highlight such pos- 
sibilities. This can be done by studying the coefficient rj 
as shown in Table I for g ^ and s = ±1. Since we have 
two possibilities for achieving bound states and scatter- 
ing, we will focus our attention first on the conditions 
giving bound states. Afterwards, only when we study 
the scattering problem is that we will take into account 
the other two conditions. 



III. SELF-ADJOINT EXTENSIONS 

For smooth functions, g G C^(M^) with g{0) = 0, we 
should have 



Hg = Hog, 



(36) 



and hence it is reasonable to interpret [48-50] the Hamil- 
tonian (31) as a self-adjoint extension of 



H 



^\{o})- 



(37) 



In order to proceed to the self-adjoint extensions of (32), 
we decompose the Hilbcrt space — L^(R^) with respect 
to the angular momentum Sj = Sjj. ® S^ip, where S^^ = 
L'^{R'^,rdr) and Sj^ = L^{S^,dip), with denoting the 



unit sphere in R^. The operator —9^ is essentially self- 
adjoint in L'^{S^,d(p) [17] and we obtain the operator 
"Ho in each angular momentum sector. Now, using the 
unitary operator V : L'^{M.'^ ,rdr) L^(IR"'", dr), given 
by {Vg){r) ~ r^l'^g[r\ the operator "Hq becomes 



YUoV' 



1 

2M 



dr^ 



+ J 



1 



(38) 



which is essentially self-adjoint for |j| > 1, while for 
|j| < 1 it admits a one-parameter family of self-adjoint 
extensions [17]. To characterize this family, we will use 
the KS [35] and the BG [36] approaches, both based in 
boundary conditions. 

In the KS approach, the boundary condition is a match 
of the logarithmic derivatives of the zero-energy solutions 
for Eq. (30) and the solutions for the problem "Ho plus 
self-adjoint extension. In the BG approach, the bound- 
ary condition is a mathematical limit allowing divergent 
solutions for the Hamiltonian (32) at isolated points, pro- 
vided they remain square integrable. 



A. KS method 

Now, the goal is to find the bound states for the Hamil- 
tonian (31). Following [35], we temporarily forget the b 
function potential and find the boundary conditions al- 
lowed for "Hq. But the self-adjoint extension provides an 
infinity of possible boundary conditions, so that it can 
not give us the true physics of the problem. Neverthe- 
less, once the physics at r = is known [3, 51, 52], it 
is possible to determine any arbitrary parameter com- 
ing from the self-adjoint extension, in such a way it is 
possible to obtain a complete description of the prob- 
lem. Since we have a singular point, we must guarantee 
that the Hamiltonian is self-adjoint in the region of mo- 
tion. Note that even if = "Hqi their domains could 
be different. The operator "Hq, with domain 'D{Ho), is 
self-adjoint if I?('HJ) = VCHo) and = 'Ho. We must 
find the deficiency subspaces, Af±, which are defined by 
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[17] 

A/"- = { v e PChJ), -hSv ^-V', 3 z- < 0} , 

with dimensions n+ and n_, respectively, which are 
caUed deficiency indices ol Ho [17]. A necessary and suffi- 
cient condition for being essentially self-adjoint is that 
n+ = n_ = 0. On the other hand, if n_|_ = n_ > 1 the 
operator Hq has an infinite number of self-adjoint exten- 
sions parametrized by the unitary operators U : Af+ — > 
AA_. 

Next we substitute the problem in Eq. (30) by the 
eigenvalue equation. 



(39) 



with fg^£ labeled by a parameter g which is related to the 
behavior of the wave function in the limit r — > rg . But for 
imposing any boundary condition (e.g. / = at r = 0) 
it is necessary to discover what boundary conditions are 
allowed to Hq. 

Now, in order to find the full domain of "Ho in 
L^( ,rdr), we have to find its deficiency subspace. To 
do this, we solve the eigenvalue equation 



(40) 



where fco G M is introduced for dimensional reasons. 
Since Hq = Hq, the only square integrable functions 
which are solutions of Eq. (40) are the modified Bessel 
functions, 



(41) 



with e = 2iMko. These functions are square integrable 
only in the range j G (—1,1), for which Hq is not self- 
adjoint. The dimension of such deficiency subspace is 
(n+,n_) = (1, 1). So we have two situations for j, i.e.. 



1< J < 0, 
< j < 1. 



(42) 



To treat both cases of Eq. (42) simultaneously, it is more 
convenient to use 



(43) 



Therefore, according to the von Neumann-Krein the- 
ory, the domain of Hq is given by 



v{hI) ^v{Ho)®^f+®^f-■ 



The self-adjoint extension approach consists, essentially, 
in extending the domain 'D{Ho) to match 2?('Ho) and 
therefore turning Ho self-adjoint. We then get 



v{no.g) = v{nl) = v{no)®N+®j^-, 



(45) 



For each g, we have a possible domain for Hq. But it will 
be the physical situation the factor that will determine 



the value of g. Thus, D{Ho^g) in _L^(R+, rdr) is given by 
the set of functions [17] 

faAr) = fu\ir) + C [i^|,|(rV^) + e^^Ky^{rVe)] , (46) 

where f\j\{r), with /|j|(ro) = f\j\{ro) = (/ = df/dr), is 
the regular wave function when we do not have Z^short('')- 
The last term in Eq. (46) gives the correct behavior for 
the wave function when r = tq. The parameter g £ [0, tt) 
represents a choice for the boundary condition. As we 
shall see below, the physics of the problem determines 
such parameter without ambiguity. In fact, g describes 
the coupling between Z^short('') and the wave function. 
Thus, it must be expressed in terms of a, the defect 
core radius rg and the effective angular momentum j. 
The next step is to find a fitting for g compatible with 
Wshoit(?')- In this sense, we write Eqs. (30) and (39) for 
f = 0, 



d^ Id f - 

—r^ H ; 7 + short 

dr^ r dr 



/o = 0. 



d^ 

dr'^ 



1 d 

r dr 



fp.o = 0, 



(47) 



(48) 



implying the zero-energy solutions /o and /p,o- Now, we 
require the continuity for the logarithmic derivative 



/o 
Jo 



fe,o 



(49) 



The left-hand side of this equation can be achieved inte- 
grating (47) from to ro. 



1 d 
r dr 



dh{r) , , 
r — -_, — I rdr = 77 



dr 



fo(r) rdr 



Mr) 



•'dr. 



(50) 



Now. considering that /o(r)/r^ does not change a lot in 
[0, ro], we find 



rdr ; 



foiro) 



/o 

So, we have 



rdr 



/o(ro) 



fo 
fo 



(51) 



(52) 



Since ro — ?> 0, the right-hand side of Eq. (49) is calculated 
using the asymptotic representation for A'|j|(x) in the 
limit X — ^ 0, given by 



-lil 



2-mr(i-|j|) 2mr(i + |j|) 



(53) 



where Cj = 7r/2 sin(7r|j|). Here a comment is in order 
upon Eq. (53). In Ref. [34], this expression has a sign 
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mistake. Here, we present the correct expression and the 
respective imphcations on the resuhs [53]. Thus, taking 
into account (46), we arrive at 



ro 



fp,0 
fp,0 



(54) 



where 



Jil 



+e 



2-mr(i-|j|) 2mr(i + |j|) 

(rVi)-l^'l (rVi)l^'l 



_2-mr(i-|j|) 2mr(i + |j|) 

Inserting (52) and (54) in (49) we obtain 



.(55) 



(56) 



which gives us the parameter g in terms of the physics of 
the problem, that is, the correct behavior of the wave 
functions when r — >■ vq, or the couphng between the 
short-ranged potential Ughovt {i^) and the wave function. 

Next we will find the bound states of the Hamiltonian 
T-Lo and using (56) the spectrum of H will be determined 
without any arbitrary parameter. Then, from Eq. (39) 
wc achieve the modified Bessel equation (k^ = —2M£, 
£ <0) 



1 d 

r dr 



h,s{r)=0, (57) 



whose general solution is given by 



(58) 



Since these solutions belong to 2?(Ho,e), it is of the form 
(46), that is, 

UAr) = fm{^) + C [A|,|(rV^) + e'^X|,|(rVi)] , (59) 

for some g selected from the physics of the problem. So, 
we substitute (58) into (59) and compute rofg,£ / fgx\r=ro 
using (53). After a straightforward calculation we have 
the relation 



feAr) 



rl^'^THM-M£)\^\ - 2I^-Ir(|j|) 



ro'''r(-|j|)(-M£)m + 2mr(b|) 



'0 



(60) 



By using Eq. (56) and solving the above equation for £, 
we find the sought energy spectrum 



2|j| + 2r7\ r(l + |j 



j2 + 21.71 + 2,7; r(i-|j|) 



(61) 



Notice that there is no arbitrary parameters in the above 
equation. 



B. BG method 

The approach used in the previews sections give us 
the energy spectrum in terms of the physics of the prob- 
lem, but is not appropriate for dealing with scattering 
problems. Furthermore, it selects the value for the pa- 
rameter g. On the other hand, the approach in [36] is 
suitable to address both bound and scattering scenarios, 
with the disadvantage of allowing arbitrary self-adjoint 
extension parameters. By comparing the results of these 
two approaches for bound states, the self- adjoint exten- 
sion parameter can be determined in terms of the physics 
of the problem. Here, all self-adjoint extensions 'Ho,Aj of 
Ho arc parametrized by the boundary condition at the 
origin. 



9o{r) = \j hm 



with 



9o{r) 



lim r'- 



'3(0- 



(62) 



(63) 



where Xj is the self-adjoint extension parameter. In [8] is 
showed that there is a relation between the self-adjoint 
extension parameter Xj and the parameter g used in the 
previous sections. The parameter g, is associated with 
the mapping of deficiency subspaces and extend the do- 
main of operator to make it self-adjoint, being a mathe- 
matical parameter. The self-adjoint extension parameter 
Xj have a physical interpretation, it represents the scat- 
tering length [54] of 'Ho,Aj [8]. For Xj = we have the 
free Hamiltonian (without the S function) with regular 
wave functions at origin and for Xj ^ the boundary 
condition in (62) permit a r"'^' singularity in the wave 
functions at origin. 



The solutions for 



Bound states 



(64) 

with fc^ = 2M£ for r 7^ 0, taking into account both cases 
in (42) simultaneously, can be written as (p ~ 2ikr) 

Uir) = A,e-«pl^lM(i + 1 + 2\j\,p) 

+i?,e-«p-l^lA/(i-|j|,l-2|j|,p), (65) 

where M(a, h, z) represents the confluent hypergeometric 
function of the first kind, and Aj , Bj are the coefficients 
of the regular and irregular solutions, respectively. By 
implementing Eq. (65) into the boundary condition (62), 
we derive the following relation between the coefficients 



Aj and Bj-. 



X,A, = [2ik)-^y^\BAl 



4(l-|j|)r^o+ 



lim r2-2bl j. (66) 



7 



In the above equation, the coefRcient of Bj diverges as 
hmr_^o+ I if lil > 1- Thus, Bj must be zero for 
\j\ > 1 , and the condition for the occurrence of a singular 
solution is \j\ < 1. So, the presence of an irregular solu- 
tion stems from the fact the operator is not self-adjoint 
for \j\ < 1, and this irregular solution is associated with 
a self-adjoint extension of the operator Ho [55, 56]. In 
other words, the self-adjoint extension essentially con- 
sists in including irregular solutions in I? ("Ho), which al- 
lows to select an appropriate boundary condition for the 
problem. 

In the present system the energy of a bound state has 
to be negative, so that fc is a pure imaginary, k ~ in. 
With the substitution k — ^ ik, we have (p' = — 2Kr) 



fi{r)=A,e-^p' 



lil 



A 



l + 2|j|,p') 
(^-|j|,l-2b|,p')- (67) 



For Eq. (67) representing a bound state, the solution 
/g (r) must vanish for r — )■ oo, i.e., it must be nor- 
malizable. By using the asymptotic representation of 
M(a, 6, z) for r — > oo, 



M(a,6;z) 



m 

r(a)' 



m 

r(b-a) 



(68) 



where T(z) is the Euler T function, the normalizability 
condition yields the relation 



A 



ni + 2\j\) ^ r(i-2|j|) 



= 0. 



Or. using the well-known relations for the T function 



(69) 



r(z-l-l) = zr(z), r(2z) = (27r)" 



^r(z)r(z- 



(70) 



we can write 



-16l' 



1^(1 



r(i 



From Eq. (66), for |j| < 1 we have 



and by using Eq. (71), the bound state energy is 



1 r(i + bl) 



1 i/iji 



A,r(i-bl) 



This coincides with Eq. (3.13) of Ref. [20] for a 
i.e., the spin-1/2 AB problem in Euclidean space, 
comparing Eq. (73) with Eq. (61), we find 



(71) 



(72) 



(73) 

= 1, 
By 



1 

a" 



1 



2|j| 



2\j\ + 2r, 



P + 2\j\+27j 



(74) 



We have thus attained a relation between the self-adjoint 
extension parameter and the physical parameters of the 



problem, j and tq. It should be mentioned that some 
relations involving the self-adjoint extension parameter 
and the S function coupling constant were previously ob- 
tained by using Green's function in Ref. [21] and the 
renormalization technique in Ref. [19], being both, how- 
ever, deprived from a clear physical interpretation. 
The (unnormalizcd) bound state wave function is given 

by 



2m 
WW) 



K\j\{-V-2M£ r)e*('"-''/2)'^, (75) 



where Ky^{z) is the modified Bessel function of the sec- 
ond kind. 



2. Scattering 

Once the bound energy problem has been examined, 
let us now analyze the AB scattering scenario. In this 
case, the boundary condition is again given by Eq. (62), 
but with the replacement Xj — )■ A|, where A| is the self- 
adjoint extension parameter for the scattering problem. 
In the scattering analysis it is more convenient to write 
the solution for Eq. (64) in terms of Bessel functions 



f£{r)^C,J{\:jlkr) + D,Yi\jlkr), 



(76) 



with Cj and Dj being constants. Upon replacing /^(r) 
in the boundary condition (62), we obtain 

Xpj^k^i^ ^ Dj[Ck-^^^-X'j{/3k^^^+C,^k-^^^ lim r^"^'^'')], 

(77) 



where 



2iJir(i + |j|)' 



C = - 



2iJ-ir(|j 



cos(^|j|)r(-|j|) k^ 

^ 7r2m ' 4(1 -|j I)- 



(78) 



As in the bound state calculation, whenever \j\ < 1, we 
have Dj ^ 0; again, this means that there arises the con- 
tribution of the irregular solution Y at the origin when 
the operator is not self-adjoint. Thus, for \j\ < 1, we 
obtain 



(79) 



and by substituting the values of j, ( and 77 into above 
expression we find 



(80) 



where 



AU^2mr(i-|j|)sin(|j|^) 



A^pb-ir(i - cos(^IjI) + 4b-ir(i + IjI) ■ 
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Since 5 is a short range potential, it follows that the 
behavior of | for r — > oo is given by [57] 



cos kr - ^ --+ 5/{k, 
TTKr V 2 4"' 



(82) 

where S^^ {k, (f>) is a scattering phase shift. The phase 
shift is a measure of the argument difference to the 
asymptotic behavior of the solution J(|m|, kr) of the ra- 
dial free equation that is regular at the origin. By us- 
ing the asymptotic behavior of J(:/, x) and Y{i>, x) for 
r — >■ oo, given by [58] 



J(v, x) 

and 

Y{,y,x)r 
respectively, we obtain 



2 / my TT 

cos X 

nx \ 2 A 



2 . / TTl' 7r^ 

sm X 

TTX \ 2 4. 



(83) 



(84) 



cos k 



-u, sm kr 

\ 2 4 



(85) 



By comparing the above expression with Eq. (82), we 
have 



cos k 



Tr\j\ TT 



a; = 



cos{kr-^-j+s'hk:cj^)]m 



with 9\s given as 



cos^Aj = dj, sin^^^i = C'j/iJ^ . 
Therefore. Eq. (82) is satisfied if 

-1/2 



1 +(.;')• 



(88) 



Now, comparing the arguments of the cosines above, the 
following phase shift is achieved: 



where 



and 



A™((/') = ^(|m|-|m-f 01), 



6'a« = arctan {fi,'] 



(89) 



(90) 



(91) 



Therefore, the scattering operator S-' . (S matrix) for the 
self-adjoint extension is 



<f>,j 



that is, 



g2iA,„(0)g2»e^ 



(92) 



4>J 



(93) 



Using Eq. (81), we have 

Ajfc2|jlr(l - |j|)e*IJ> -I- 4iJ'ir(i + |j|) 



A^fc2|ilr(l - |j|)e-'IJ> + 4lilr(l -I- |j|) 



(94) 

In accordance with the general theory of scattering, the 
poles of the S matrix in the upper half of the complex 
plane [59] determine the positions of the bound states in 
the energy scale, Eq. (73). These poles occurs when the 
denominator of Eq. (94) is equal to zero and with the 
replacement k — > iK. So, we have 

A^'(iK)2|jlr(l - |j|)e-'IJ'l" + 4l^lr(l + |j|) = 0. (95) 

Solving this equation for £, we found 

i/IjI 



1 r(i + |j|) 
K r(i-|j|) 



(96) 



Thus, by comparing the expression above with Eq. (73), 
we find = Aj, with A^- given by Eq. (74), and the self- 
adjoint extension parameter for the scattering scenario 
being the same one as that for the bound state problem. 
This is a very interesting result first discussed in [34]. 
Thus, we also obtain the phase shift and the scattering 
(87) matrix in terms of physics of the problem. If AJ = 0, we 
achieve the corresponding result for the pure AB prob- 
lem with Dirichlet boundary condition; in this case, we 
recover the expression for the scattering matrix found in 



Ref. [60], S, 



'. If we make Aj = oo, we get 



„2iA„(0)-|-2i7r|j| 



(97) 



The scattering amplitude f-^ ^{k, 6) can be obtained using 
the standard methods of scattering theory, namely [fk = 




2iA,- 



1 e 



imO 



1 










1 
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IV. AHARONOV-BOHM PROBLEM PLUS A 
TWO-DIMENSIONAL HARMONIC OSCILLATOR 

In this section, we address the AB problem plus a two- 
dimensional harmonic oscillator (HO) located at the ori- 
gin of a polar coordinate system. The potential of the 
harmonic oscillator in two-dimension space is given by 



Vho - + ]^Mujy 



In polar coordinates (r, tp) it can be written as 



Vho = ^MwV^, 



(98) 



(99) 



where we are considering an isotropic HO, lOx = = t^- 
Let us now include the oscillator potential (99) into 
the AB Hamiltonian (25), which leads to the following 
eigenvalue equation 

Hho^ = «5t$, (100) 

where 

Hho=Hnr + Vho. (101) 
By writing the solutions of Eq. (100) in the form 



we obtain the radial equation 

n'<i>e,{r)^£'<j>e'{r), 
where H' = Wq + Z^short, 



1 

2M 



1 d 
r dr 



2 2 

7 r 



(102) 



(103) 



(104) 



and 7^ = M^uj^. 

In order to have a more detailed analysis of this prob- 
lem, we will first examine the motion of the particle con- 
sidering two cases (i) excluding the r = region and (ii) 
including the r = region afterwards. At the end, we 
compare with some results in the literature. 



3. Solution of the problem excluding r — region 

In this case, the Hamiltonian (101) does not contain 
^short- By direct solving 



H^,0£,(r) (r). 



(105) 



we obtain 



1+3 



i£'(r) = Aj7— r^e~27'- M{d,l+ j.-fr^) 



1+1 4 _l-v^2 



where d = (l+j)l2 - M£' /2-i, U{a,b,z) is the conflu- 
ent hypergeometric function of the second kind, and Aj, 
Bj arc constants. However, as only M is regular at ori- 
gin, it should be imposed Bj = 0. Moreover, if d is 

or a negative integer, and 1 + |j| is a negative and 
non null integer, the series terminates and the hyperge- 
ometric function becomes a polynomial of degree n [58] . 
This condition guarantees that the hypergeometric func- 
tion is regular at origin, which is essential for the treat- 
ment of the physical system since the region of interest 
is that around the charge filament. Therefore, the series 
in (106) must converge if we consider that d = —n and 

1 + IjI where n,£ e Z*, with Z* denoting the set of 
the nonnegative integers. This condition also guarantees 
the normalizability of the wave function. So, using this 
condition, we obtain the discrete values for the energy 
whose expression is given by 



£' = {2n + l + \3\)u 



(107) 



where n e Z and 1 + |j| 7^ — ^, ^ G Z. The (unnormal- 
ized) energy wave function is given by 



$(r, if) 



I, -I 

7 2 r"'e 



i7r^gi(m-s/2)ip 

xA/(-7i,l + |j|,7r2), (108) 



where Cj is a normalization constant. Notice that in 
Eq. (107) |j| can assume any non integer number larger 
than —1. However, we will see that this condition is no 
longer true when we include the term lAshon- Next to 
study the motion of the particle in all space, including 
the r = region, we invoke the self-adjoint extension of 
the operator Hg. 



^. Solution including the r = Q region 

In this case, the dynamics includes the Wshort term. So, 
let us follow the same procedure as in Sec. Ill A to find 
the bound states for the Hamiltonian T-L' . Like before we 
need to find all the self-adjoint extension for the operator 
Hq. The relevant eigenvalue equation is 



H;</'^,£'(r) =f>^,£'(r). 



(109) 



with (pi,^ labeled by a parameter -d. The solution to this 
equation are given in (106). However, the only square 
integrable functions is U{d, 1 +i, jr'^). Then, this implies 
that Aj = in Eq. (106). Therefore, we get 



^'^'^ U{d,l+j,-fr^). 



(110) 



In order to guarantee that (l>£'{r) G i^(R, rdr), it is ad- 
visable to study their behavior as r — > 0, which implies 
analyzing the possible self-adjoint extensions. 

Now, to construct the self-adjoint extensions, let us 
consider the eigenvalue equation 



+Bj7— r^e-^-^"^ C/(d,l+j,7r^), (106) 



±iko(j>±{r). 



(Ill) 
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Since Hq = H'q and, from (110), the solution to this 
equation is given by 



(112) 



where d± = (1 + =F iko/A-y. 

Let us now consider the asymptotic behavior of 
;7(d±,l+j,7r2) asr^O [61] 



C/(d±,l+j,7r2) 



r(rf±) r(d±-j) 



(113) 



Working with this expression, let us find under which 
condition the term. 



|(/)i^±(r)|Vdr, 



(114) 



has a finite contribution near origin region. Taking Eqs. 
(112) and (113) into account, we have 



lim |0i.±(r)| r 

r— >0 



[A 



1+2] 



2Jl 



(115) 



where ^land A2 are constants. Studying Eq. (115), we 
see that <j)±{r) is square-integrable only for j e (—1,1). 
In this case, since is expanded by 4>+{r) only, we 
have that its dimension is rt+ = 1. The same applies to 
N- and 4>-{r), resulting in n_ = 1 . Then, H'q possesses 
self-adjoint extensions parametrized by a unitary matrix, 
U(l) = e'^, with I? e [0,27r). Therefore, the domain of 



Til is given by 



(116) 



So, to extend the domain T>{Hq) to matchI?(H[/) and 
therefore to turn Hq self-adjoint, we get 

ViHi^) = V{H'^) = 2?(H[,) © A/-+ ® AA_, (117) 

we mean that, for each t?, we have a possible domain for 
'D{'H'q ^), but it will be the physical situation which will 
determine the value of -d. The Hilbert space (117), for 
both cases of Eq. (42), contains vectors of the form 



(b^Ar) = 'Pljlir) + C rl^le-— C/(d+, 1 + \j\,^r') 



+e 



;7(d_,l + |j|,7r')]}, 



(118) 



where C is an arbitrary complex number, (/)|j|(0) = 
4>\j\{Q) = and 4'\j\{r) g L'^{M+ ,rdr). For a range of -d, 
the behavior of the wave functions (118) was addressed in 
[62] . However, as we will see below, if we fix the physics of 
the problem at r = Tq, there is no need for such analysis 
because the value of '9 is automatically selected. 

In order to find the spectrum of "Hg we consider the 
limit r — >■ of <^^,£'(r), i.e., we substitute (110) in the 
left-hand side of (118) and, using (113), after equating 
the coefficients of the same power in r leads to 



A 

m 



• 1 e"^ 

.f(d7) ^f(d^ 



(119) 



and 



A 



1 



AS 



T{d^\:i\) "nd+ 

whose quotient leads to 



iji) nd- - iji 



:), (120) 



1 



T{d) 



(121) 



r(d+-b|) ' r(d_ 



The left-hand side of this equation is a function of the 
energy £' while its right-hand side is a constant (even 
though it depends on the extension parameter i? which is 
fixed by the physics of the problem) . Then we have the 
equation 



r ( _ M£L 

V 2 27 

r ( _ M£L\ 

V 2 27 ; 



const. 



(122) 



Therefore, we have achieved the energy levels with an 
arbitrary parameter i?, and we get inequivalent quantiza- 
tions to different values of it [63] . Each physical problem 
selects a specific value of 

Next we replace problem defined by the Hamiltonian 
%' = Hq +Wshort, by H'q plus self-adjoint extensions. 
So, we consider the static solutions {£' = 0) 4>o{r) and 
(jjsfiif) for the equations 



(123) 



(124) 



Now, to find out the value of it is required the conti- 
nuity of the logarithmic derivative 



Mr) 



Mr) 



4>Sfiir) 



4>s,o{r) 



(125) 



The left-hand side of Eq. (125) is obtained by integration 
of Eq. (123) from to tq. Since the integration of the 
harmonic term 



7- / r 
/o 



ro pro 

3+2|j|</)o(r)drW720o(ro) / r^+^^'lrfr 







0, 

(126) 

as ro — >■ 0, the result is the same as in Eq. (52). Then, 
we arrive at 



Mr) 



Mr) 



(127) 



The right-hand side of Eq. (123) is carried out using 
the asymptotic behavior of U{d, 1 -I- |j|, 7r^), given in Eq. 
(113), in Eq. (118), and it take the form 



4's,o{r) 
(I^Sfiir) 



%{r) 



(128) 
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where 



-\3\ 



n-\j\)r 



+e 



Tid+) r(d+-|j|) 

n\j\)hr)-y^ , r(-|j|)rl^-| 



(129) 

Replacing Eqs. (127) and (128) in expression (125), we 
obtain 



(130) 



With this relation, we select an approximated value for 
parameter t9 in terms of the physics of the problem. 

The solutions for (105) are given by (110) and since 
this function belongs to 'D{Hq^), it is the form (118), 
that is. 



U{d+,l + \j\nr^) 



+e'''C/(d_,l + b1,7r')]} 



(131) 



for some selected form the physics of the problem at r = 
ro. So, using (110) in (131) we have (6 T{d) /T{d-\j\)) 



To 



r^i%(-|j|)7i^"ie-r(b|) 



-|j|)7i^-ie + r(|j|) 



(132) 



Using the Eq. (130) in the above equation we obtain 



j2-2b-| + 2rA r(l + |j|) 



r(d - ^b\rf^\ \.f + 2\j\ + 2r,J T{1 - ' 

(133) 

The Eq. (133) is too complicated to evaluate the bound 
state energy explicitly, but its limiting features are inter- 
esting. If we take limit ro — > in this expression, the 
bound state energy are determined by the poles of the 
gamma functions, i.e.. 



-1< J <0, £' 
< J < 1, £' 



(2n 
(2n 



£' = {2n + l±\j\)o 



(134) 



(135) 



where n <S Z*. The + (— ) sign refers to solutions which 
are regular (irregular) at the origin. This result coincide 
with the Eq. (1) of Ref. [64], for the special case of 
a ~ 1. Another interesting case is that of vanishing 
harmonic oscillator potential. This is achieved using the 
asymptotic behavior of the ratio of gamma functions for 
7^0 [65], 



r(^-ff) 

r { LM _ ml) 

V 2 27 ; 



27 ; 



(136) 



which holds for f < and this condition is necessary for 
the pure AB system to have a bound state. Using this 
limit in the Eq. (133), one finds 



-MS 



1 



.72-2|j| + 2r;\ r(l + |j|) 



^bV^I^I \f + 2\j\ + 2rjj T{l-\j\)- 

(137) 

Then, by solving Eq. (137) for £, we obtain 



Mr2 



r(l + |j|) f.f-2\j\ + 2y 



r(l-|j|) \ f + 2\j\ + 2T^ 



iI/IjI 



(138) 



which is the result obtained in Eq. (61). Thus, in the 
limit of vanishing harmonic oscillator, we recover the 
pure AB problem. 

Now we have to remark that this result contains a sub- 
tlety that must be interpreted as follows: the presence of 
the singularity in the problem establish the range of j 
as given by |j| < 1. If we ignore the singularity and 
impose that the wave function is regular at the origin 
((/)(r) = 0(r) = 0), we achieve the same spectrum of Eq. 
(135), but with j assuming any non integer value greater 
than —1 [66-68]. In this sense the self-adjoint exten- 
sion prevents us from obtaining a spectrum incompatible 
with the singular nature of the Hamiltonian when we 
have Z^short [69, 70]. We have to take into account that 
the true boundary condition is that the wave function 
must be square-integrable through all space, regardless 
it is irregular or regular at the origin [35, 70]. 



5. Determination of self-adjoint extension parameter 

For our intent, it is more convenient to write the solu- 
tions for (103) for r 7^ 0, taking into account both cases 
in (42) simultaneously, as 



(P£'[r) = Ajj 2 e 



>l^lAf(d,l + ljl,7r2) 

^r-l^lA/(d-lj|,l-ljl,7r2), 

(139) 



where Aj , Bj are the coefficients of the regular and sin- 
gular solutions, respectively. By implementing Eq. (139) 
into the boundary condition (62), we derive the following 



relation between the coefficients Aj and Bj : 



A'A,7l^l=S, 1- 



£' 



4(1 -IjI) 



lim r 



(140) 



where A^- is the self-adjoint extension parameter for the 
AB plus HO system. In the above equation, the coeffi- 
cient of Bj diverges as limr_i.o+ r^"^!-'!, if \j\ > 1. Thus, 
Bj must be zero for \j\ > 1, and the condition for the oc- 
currence of a singular solution is [j] < 1. So, the presence 
of an irregular solution stems from the fact the operator 
is not self-adjoint for \j\ < 1, recasting the condition 
of non-self-adjointness of the previews sections, and this 
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irregular solution is associated with a self-adjoint exten- 
sion of the operator H'q [55, 56]. In other words, the 
self-adjoint extension essentially consists in including ir- 
regular solutions in 2?('Ho) to match 2?('Hq^), which al- 
lows us to select an appropriate boundary condition for 
the problem. 

In order to Eq. (139) to be a bound state, (j>j^£' (r) must 
vanish for large values of r, i.e., must be normalizable 
at large r. By using the asymptotic representation of 
M(a, b, z) for z — > cx), 

MiaAz)-^^fz-^^^^i-z)-^ (141) 

the normalizability condition yields the relation 



r(i + l,l)r(^ 



Me' \ 

27 ) 



r(i-|j|)r(i^-M^ 

I 2 27 



A. 



(142) 



From Eq. (140), for |j| < 1 we have = A^^I-^Uj and 
by using Eq. (142), the bound state energy is implicitly 
determined by the equation 



r 



1+1 J I 



MS' 

27 



1 r(i + |j|) 



i-IjI 



27 



a;.7I^i r(i 



1.71 



(143) 



This results coincides with Eq. (53) of Ref. [21] for a = 1, 
as it should be. By comparing Eq. (143) with Eq. (133), 
we find 



1 ^ 2 / j2-2|i| + 277 



(144) 



We have thus attained a relation between the self-adjoint 
extension parameter and the physical parameters of the 
problem. 

The limiting features of Eq. (143) are interesting. For 
Aj — > or CX3, from the poles of gamma function, we have 



A' 



0, £' = (2n+l + |j|)c^, 
oo, f = (271 + 1 - \i\)uj. 



(145) 



These bound states energies coincide with those regular 
and irregular solutions given in Eq. (135) of the previous 
section. From relation (144) the regular wave function, 
when A^ = and the spin is absence, it is associated with 
< J < 1, and the irregular wave function, when A^ = 
oo, is associated with —1 < j < 0. Another interesting 
limit is that of the vanishing oscillator potential. So, 
using the limit (115) in (73), we have 



1 r(i + |j|) 
A,7i^ir(i-|j|)- 



(146) 



with \j the self-adjoint extension for the pure AB system. 
From this equation we have the spectrum for the pure AB 
system in terms of the self-adjoint extension parameter. 



£• = - 



M 



1 r(i + |j|) 



A,r(i-|j|) 



(147) 



This result also coincides with Eq. (3.13) of Ref. [20] for 
a = 1. By comparing this equation with (138) we arrive 
at 



1^_J_ / j^~2|j| + 2?/ 
A,- rJ}o\\p + 2\j\ + 2r^ 



(148) 



for the relation of the self-adjoint extension parameter 
and the physics of the problem for pure AB system. 
Then, the relation between the self-adjoint extension pa- 
rameter and the physics of the problem for the pure AB 
has the same mathematical structure as for the AB plus 
HO. However, we must observe that the self-adjoint ex- 
tension parameter is negative for the pure AB, confirm- 
ing the restriction of negative values of the self-adjoint 
extension parameter made in [21], in such way we have 
an attractive 5 function. It is a necessary condition to 
have a bound state in the pure AB system. 



V. CONCLUSIONS 

We have presented the detailed calculations of Ref. 
[34], in which was proposed a general regularization 
method to address a system endowed with a singular 
Hamiltonian (due to localized fields sources or quantum 
confinement). Using the KS approach, the bound states 
were determined in terms of the physics of the problem, 
in a very consistent way and without any arbitrary pa- 
rameter. In the sequel, we employed the BG approach. 
By comparing the results of these approaches, we have 
determined an expression for the self-adjoint extension 
parameter for the bound state problem, which coincides 
with the one for scattering problem. We have thus ob- 
tained the S matrix in terms of the physics of the problem 
as well. In this point, we remark that the important re- 
sults of Refs. [3, 20, 71] are giving in terms of the an 
arbitrary self-adjoint extension parameter. In our work 
this parameter was determined in terms of the physics 
of the problem. The outcomes obtained by Park are a 
particular case of our results for a fixed value of the self- 
adjoint extension parameter. To our knowledge, it was 
not known in the literature an expression for the bound 
state energies for the AB with a defined self-adjoint ex- 
tension parameter. In Ref. [34] this expression was pre- 
sented by the first time, whose details are derived here. 

To illustrate the applicability of our approach to other 
physical systems, we deal with the spin- 1/2 AB problem 
plus a two-dimensional harmonic oscillator. Two cases 
were considered: (i) without and (ii) with the inclusion of 
short range potential Z^short in the nonrelativistic Hamil- 
tonian. Even though we have obtained an equivalent 
mathematical expression for both cases, has been shown 
that, in (i) j can assume any noninterger value greater 
than —1 while in (ii) it is in the range j E (—1, 1). In the 
first case, it is reasonable to impose that the wave func- 
tion vanish at the origin. However, this condition does 
not give a correct description of the problem in the r = 
region. Therefore, the energy spectrum obtained in the 
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second case is that physically acceptable. The presence 
of the singularity establishes that the effective angular 
momentum must obey the condition j G (—1,1) and im- 
plies that irregular solutions must be taken into account 
in this range. The only situation in which we can neglect 
the short range potential is the one in which one looks 
only for topological phases. 

A natural extension of this work is the inclusion of the 
Coulomb potential [72], which naturally appears in two- 
dimensional systems, such as graphene [73] and anyons 
systems [74, 75]. 
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